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1. Introduction 

The simplest lattice discretization of the S\J{Nc) Yang-Mills theory is the well-known Wilson 
plaquette action [1]. Different discretization of the lattice action will not change the physics of 
the continuum limit realised in the neighborhood of the weakly coupled ultravioled fixed point. 
However, far from this continuum limit, different discretizations can lead to the appearence of 
second order phase transition points that can mimic a continuous infrared fixed point for the theory 
defined by the naive lattice dicretization. In fact, although a continuum theory can be defined at 
any of those points, in principle this theory is not related to the ultraviolet gaussian fixed point. 
One possible extension of the Wilson action includes plaquette terms in a representation of the 
gauge group other than the fundamental. For example the following action includes a term in the 
adjoint representation 

S = iSfund L f 1 - ^Re Trf {UMi))) +i8adj L f 1 " 77237^^ ^r^ {UMi))) . d-D 

!,Ai>v V ■'V / ,;^>v V ■'V ^ / 

where Nc is the number of colours and U^y{i) the plaquette in the (jLt, v)-plane from point /. The 
sum over all the points / is done over the four-dimensional hypercubic lattice L^. Tip and Tr^ are, 
respectively, the trace defined in the fundamental and in the adjoint representation of the SUCA'^c) 
gauge group. They are related by Tr a{U) = \Tr f{U)\'^ — 1. 

This fundamental-adjoint plaquette action has been used in the pioneering work of Ref. [2] and in 
several more recent studies [3] which extensively investigated the structure of the phase diagram 
for Nc = 2. Our interest in this model comes from the recent studies of the conformal window using 
lattice field theory techniques. The SU(2) gauge theory with 2 adjoint fermions has been shown 
to have an infrared conformal fixed point by looking at the scaling properties of the mesonic and 
gluonic spectrum [4]. In principle, the same features could be reproduced around a second-order 
phase transition point appearing as a lattice artefact due to the chosen lattice discretization. The 
fundamental-adjoint lattice model of Eq. (1.1) can be seen as the leading contribution to the action 
with adjoint fermions in the heavy bare quark mass limit. Hence, if the end-point of the first order 
phase transition in this model turns out to be a lattice-induced second order phase transition point, 
one needs to investigate how the results of Ref. [4] would be affected by it. In the following, we 
investigate carefully the phase diagram and the spectrum of the lattice model in the vicinity of the 
end-point to check whether the infrared physics resembles the one studied in Ref. [4]. A detailed 
description of our study will be the object of a forthcoming publication [5]. 

2. Phase diagram 

In the two-dimensional plane of the coupling constants (j8fund>j8adj)> the theory presents several 
regions (cfr. Ref. [2] for a qualitative picture). On the fundamental axis jSadj = the SU(2) gauge 
theory with standard Wilson plaquette action is recovered and this is known to have a crossover re- 
gion at jSfund ~ 2.30. When the adjoint coupling is turned on and the second term of Eq. (1.1) starts 
becoming important, the system develops a first order bulk transition which becomes stronger as 
j3adj increases. We have monitored the location of this bulk transition line by studying the expecta- 
tion value of the fundamental and adjoint plaquettes, and we have also computed the corresponding 
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normalized susceptibilities. An example of the hysteresis cycle characteristic of the bulk phase 
transition at large jSadj = 1.50 is shown in Fig. l(left): on a hypercubic symmetric lattice of size 
8^ we clearly distinguish two sepai^ate branches for the fundamental plaquette as jSfund is changed 
starting from a random (hot) or unit (cold) gauge configuration. When /3adj is decreased, we note 
that larger volumes ai^e necessary in order to con^ectly identify the presence of the hysteresis loop. 
For example, at j8adj = 1.275 a lattice 12^ is not large enough for the system to develop the two 
vacua of the first order transition, and this is shown in Fig. 1 (right). 
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Figure 1: Hysteresis cycle of the fundamental plaquette for different values of the adjoint coupling jSadj- At 
large j3adj (left), the separation between the hysteresis branches is visible on relatively small volumes. At 
smaller )3adj (right) bigger volumes are necessary to clearly distinguish the first order nature of the transition. 
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Figure 2: Fundamental plaquette difference between hot and cold start runs at the centre of the hysteresis 
cycle. Also shown is the approximate position of the critical ^adj value at which this difference is expected 
to vanish. A consistent result is found using the adjoint plaquette. 

By further decreasing jSadj, the separation between the lower and the upper branch of the hys- 
teresis shows a clear trend suggesting that it should vanish at approximately jSadj ^ 1.25. An esti- 
mate of this separation is given by the difference of the plaquette in the two vacua at the center of 
the hysteresis loop: 

Apfund = (Plaq^ , 1 > - (Plaqf ,2 > , (2-1) 
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where the subscripts 1 and 2 refer to the disctinct vacua, and a similar definition holds for A/^adj- 
This Apfund is plotted in Fig. 2, where we always used the smallest volume where the first order 
nature of the transition was manifest. The point at jSadj = 1-25 required a very large 40^ lattice for 
which we cuiTcntly do not have very good control over the systematic and statistical uncertainties 
of the simulation. 

In the region below the approximate location of the end-point, we have checked that the transi- 
tion becomes a crossover, signalled by the lack of scaling with the volume in the fundamental and 
adjoint plaquette susceptibilities. The height and the location of peak of the susceptibility is con- 
sistent across the different volumes. The location of the peak separates a strong coupling region at 
small jSfund from a region closer to the weak coupling limit (jSfund — ^ °°)- 
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Figure 3: The location of the bulk phase transition is delimited by black lines, representing the extension of 
the hysteresis cycle. The approximate location of the bulk transition end-point is shown by the arrow. The 
coloured symbols show the points where the spectrum is investigated. 

3. Spectrum measurements 

Let us first recall here that we do not want to precisely pin down the end-point location, but 
rather to identify its neighbourhood, where the spectrum of the theory should be investigated. Our 
aim is to compare the scaling properties of the spectrum when the bulk transition end-point is ap- 
proached in a controlled manner, with the ones of the model with 2 adjoint fermions. This will 
help us clarify the still controversial nature of this end-point. In Fig. 3 a summary of our results 
concerning the bulk phase transition line and its end-point is shown. In addition, we indicate the 
points where we performed a detailed investigation of the spectrum, as described in the following. 
Our simulations with the fundamental-adjoint action ai^e caiTied over using a modified Metropolis 
algorithm [6] which helps us cope with the increasing autocorrelation times due to critical slowing 
down when simulating closer to /3adj ~ 1-25. This allows us to obtain a large statistics of inde- 
pendent gauge configurations even on lai^ge lattices when the critical slowing down starts affecting 
the simulations. In particular, we measure our observables on Monte Carlo histories of ^(10000) 
configurations, each separated by i^(lOO) modified-Metropolis updates of the SU(2) link matrices 
(jSadj values closer to the end-point have a larger number of intermediate updates between mea- 
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surements to reduce autocorrelations in our ensembles). In the following, we show results at four 
different values of the adjoint coupling jSadj = 1.00, 1.05, 1.10, 1.16 and spanning a large range 
of j3fund such that both regions around the crossover are monitored. A sequence of five different 
volumes is simulated for each jSadj: 6^ x 12, 10^ x 20, 16^ x 20, 24^ x 32 and 32^ x 32, where the 
longer temporal extent is used to better identify effective mass plateaux. 

We employ the variational procedure detailed in Ref. [7] to extract the ground state mass and a few 
excitations of the spectrum in the following channels: 



String tension: oy/a is the lightest dynamical scale in a pure gauge theory and it is used to 
set the overall scale. We extract the string tension from correlators of long spatial Polyakov 
loops of length L = aNg. The asymptotic large-time behaviour of these correlators is gov- 
erned by the lightest torelon state whose mass ami can be used to obtain the string tension 
according to 

ami{N,) = a^oN, ^ . (3. 1) 



3A^, \%N^a^o 

The validity of the above equation is checked a posteriori by comparing the extracted string 
tension with the one obtained using only the leading term —n/{3Ns). Significant finite-size 
systematics are absent when L^/a > 3, which we satisfied in our simulations using large 
spatial volumes for the smallest values of a^/a. 



Scalar glueball mass: amQ++ is the lightest glueball mass in the spectrum. Correlators of 
smeared spatial Wilson loops in the scalar representation of the cubic symmetry group are 
measured on each configuration. 

Tensor glueball mass: am2++ is the second lightest glueball and its mass is monitored to 
check whether its behaviour is the same as the scalar one. Having quantum numbers different 
from the vacuum, its scaling properties could be different in principle. 
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Figure 4: Measured spectrum on different volumes at jSadj = 1 -00. (left) The square root of the string tension 
in units of the lattice spacing as jSfund increases towards the weak coupling region, (right) The mass of the 
lightest scalar glueball in units of the lattice spacing. Very good control over finite-size effects is obtained 
for both these observables. 
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In Fig. 4 we show the string tension and scalar mass at jSadj = 1 -00 and for a range of different 
Aund and spatial volumes L^. The string tension decreases monotonically when approaching the 
weak-coupling limit at large j3fund> whereas the scalar glueball mass develops a short plateaux in 
the crossover region before decreasing again. In both cases we have a good control over finite- 
size effects, with masses matching on at least two subsequent volumes for each point. The largest 
finite-size effects are seen towards the weak coupling, where the string tension becomes small. Not 
shown in the plots is the behaviour of the tensor glueball which resembles the string tension one, 
though with somewhat larger finite-volume systematics. 

Given the large range of lattice volume used, we are able to reliably estimate the infinite volume 
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Figure 5: (left) The infinite volume limit of a^/a is plotted for jSadj closer and closer to the end-point 
location, (right) The scalar glueball mass in the infinite volume limit for all available j3adj- The very different 
behaviour of the two observables is clear when /3adj is increased. Also shown are vertical dashed lines 
highlighting the location of the crossover region. 



Umit of ay/o and amQ++ for all the jSadj values studied. However, for some of these values we can 
not extrapolate am2++ at L = oo with our current data. The comparison of the extracted infinite- 
volume spectrum between different jSadj is shown in Fig. 5. To study the scaling of the observables 
along a trajectory when approaching the bulk transition end-point, we choose to follow the line 
in the phase diagram given by the peak of the fundamental plaquette susceptibility. The location 
of such peak at the four jSadj values investigated is highlighted by vertical dashed lines in Fig. 5. 
On those points, we note that the string tension remains constant when approaching the end-point 
(larger j3adj), whereas the scalar glueball mass slightly decreases: this suggest a non-constant ratio 
mo++ / y/o. A summary of this result is shown in Fig. 6. 



4. Conclusions 



In this work we have studied a SU(2) pure gauge theory with a modified lattice plaquette 
action. We added a coupling to plaquettes in the adjoint representation of the gauge group. This 
lattice system is known to have a bulk phase transition with an end-point relatively close to the 
fundamental coupling axis. The nature of this end-point is still controversial and we focused 
more on the region close to it, but fai^ from the bulk phase transition. Thanks to our improved 
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gluonic spectroscopic technique [7], we measured the string tension, the scalar glueball mass and 
the tensor one, aiming at studying their scaling properties when the end-point is approached. This 
is the first study of the gluonic spectrum in this model. Therefore we carefully checked for finite- 
size systematics and tried to reduce autocoiTclation effects on our observables. The extrapolated 
infinite-volume spectrum shows a non-constant mo++ / ^/G ratio when approaching the end-point 
in a controlled manner. This seems in contrast with the infrared dynamics of the SU(2) theory 
with 2 adjoint fermions, where such a ratio is driven by a conformal fixed point and is consistent 
with the continuum SU(2) Yang-Mills value mo++ / \/g ~ 3.7. The results presented here should 
be considered as preliminary and might be still too far from the basin of attraction of the end-point. 
More detailed analysis and discussions will be presented in a forthcoming publication [5]. 
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Figure 6: (left) The scalar glueball mass and the square root of the string tension at the values (/3fund,j3adj) 
defining the plaquette susceptibility's peak when approaching the end-point from below, (right) The ratio 
"'°t^ on the same points. We highlighted the approximate location of the bulk transition end-point. 
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